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ABSTRACT: We discuss the initial conditions problem for inflation driven by the vacuum 
energy of a plateau potential, and in particular the Starobinsky inflation. We show that 
the supergravity embedding of the R + R 2 theory naturally decreases the size of the acausal 
homogeneity, required for the low-scale inflation to occur, thanks to the presence of the dy¬ 
namical pure supergravitational “auxiliary” fields. We examine the evolution of the R + R 2 
fields within a FLRW Universe. We also find a dependence of the initial conditions problem 
on the background spatial curvature. 


ArXiv ePrint: 1502.01246 



Contents 


1 Introduction 1 

1.1 The initial conditions problem 2 

1.2 Outline of the objectives and results 5 

2 Old-minimal R + R 2 supergravity: the U sugra .R 2 8 

2.1 Evolution 10 

2.1.1 Dynamics of the expansion 13 

2.1.2 Selfreproduction 14 

2.2 Initial conditions 15 

2.2.1 R + R 2 15 

2.2.2 R + R 2 supergravity 16 

3 New-minimal R + R 2 supergravity: the C nsugrELR 2 17 

3.1 Evolution 18 

3.2 Initial conditions 20 

4 The curvature term 21 

4.1 Closed Universe 23 

4.1.1 R + R 2 23 

4.1.2 R + R 2 supergravity 24 

4.2 Open Universe 25 

5 Conclusions 27 


1 Introduction 

The hot Big Bang model provides a reliable and tested description for the cosmic evolution 
from at least as early as the time for the synthesis of the light elements, about one second 
after the Big Bang. Nevertheless, problems associated with the adiabatic expansion of the 
hot early FLRW Universe, such as the entropy, the flatness and the monopole, could not be 
explained until the postulation of the inflationary cosmology in the beginning of 1980s. [1]. 
Cosmological inflation is a phase of accelerated expansion assumed to have taken place in the 
most primal times, once the universe emerged from the Planck epoch, see e.g. [2, 3]. It has 
attracted the interest of the cosmologists for it is a theory with predictive power and can be 
naturally implemented in microscopic models [4], without the need for rather special initial 
conditions. Inflation, indeed, is the most efficient mechanism that magnifies, homogenizes 
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and isotropizes the universe and on top of that it successfully describes the CMB temperature 
fluctuations according to several data sets; the Planck 2013 and 2015 results [5-8] is the 
latest example. Inflation can be implemented simply by a theory whose matter content acts 
as vacuum energy, or equivalently by a quasi cosmological constant term, and magnitude 
comparable to the GUT scale. 

PLANCK CMB data [6, 8], although strongly support the basic picture of the inflationary 
theory may question the generality of the cosmological phase because special initial conditions 
seem to be required [9-13]. The absence of primordial tensor modes 1 and the spectral index 
values for the scalar perturbations favour the plateau-like potentials. These potentials are 
characterized by relatively low energy densities and are not capable to drive an inflationary 
phase right after the Planck era implying that our Universe started with a decelerating phase 
(& < 0) instead of an accelerating one (a > 0). It is known that low energy scale inflation 
renders the generality of the initial conditions subject to speculation, see e.g. Goldwirth and 
Piran [15] for a classical review on this topic. Motivated by these observations, in this work, 
we examine the initial conditions required for the R + R 2 gravity and supergravity models of 
inflation. 

1.1 The initial conditions problem 

In the Friedman-Lemaitre-Robertson-Walker (FLRW) model the Universe has a finite age and 
the cosmological scale factor a{t) grows slower than t, facts which lead to the notion of the 
horizons: the particle and the event horizons that grow linearly with time. Hence, regions in 
the universe have past histories that do not interact. The present observable universe consists 
of about 10' 5 causally disconnected regions at the epoch of recombination. However, the CMB 
is uniform to less than one part in 10 . It has been shown that the class of all initial conditions 
for which the universe at late times behaves as an FLRW Universe is of measure zero, see e.g. 
Collins and Hawking [16]. It is unlikely that the universe began in a chaotic state and has 
reached the CMB homogeneous state in the course of an adiabatic evolution. This reasoning 
motivated the introduction and establishment of the inflationary theory. 

According to the latest PLANCK data 2 , the inflationary models fully consistent with the 
data are the plateau-like potentials with a representative example the Starobinsky model [17]. 
It is an f{R) = R/2 + R 2 /(12m 2 ) gravity theory which in the dual picture yields a potential 
for the scalaron tp that reads 

Vi&iv) = ^inf ^1 — e \R ,p ! Mp \ . (1.1) 

The Linf = (3/4)m 2 Af|, is the characteristic upper bound of the inflationary energy density 
for the Starobinsky model, U(nf ~ 1CU 10 MJ, <C Mp. This model, although it originally 
accounts for one of the first attempts to describe the evolution of the universe in its earliest 

1 Here, we consider that the signal of BICEP2 CMB experiment is contaminated by foreground dust [14]. 

2 The discussion and the results of the paper were based on Planck 2013 data. This version (v2) also cites 

the Planck 2015 data which further support the plateau inflationary potentials. 
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Figure 1: In the chaotic models (left panel), inflation starts from Planck densities. In the Starobin- 
sky I? 2 -model (right panel) inflation cannot start at the Planck time because the inflationary sector 
of the potential is bounded from above, Vp? < Finf •C Mp. 


moments, it differs from the dominant pre-PLANCK inflationary picture. For decades the 
standard paradigm for inflation has been that our Universe emerged from the quantum gravity 
era, and has been magnified to cosmological scales thanks to the prevailing presence of the 
potential energy of a scalar field in the energy-momentum tensor. The common illustration of 
this paradigm has been the quadratic large field model, V(4>) = m 2 (j) 2 , which is at the edge of 
the 95% CL contours allowed by Planck WP-sTiigh-£ CMB data [6] (and ruled out at over 
99% according to Planck 2015 data sets [7]). 

After PLANCK, inflation appears to start at low energies within a pre-existing homo¬ 
geneous initial patch. Hence, from one point of view, the plateau potentials bring back the 
problematic requirement of initial acausal homogeneity a fact that renders inflation nongeneric 
and reduces some of its appealing power to free cosmology from the need for specific initial 
conditions. 

It is very motivated the inflationary phase to have been initiated close to the Planck 
energy scales for it assures the natural creation of our observable Universe without rather 
special initial conditions. Indeed, even a fundamentally small initial patch of Planck length 
radius Ip when dominated by the potential energy of the inflaton field, ^<f 2 + ^(Vc/>/a) 2 < 
V(<fi) ~ ptot ~ Mp, starts expanding in an accelerating manner. The essential implication of 
this accelerated expansion is the presence of a nearly constant event horizon distance whose 
size is also ~ Ip, that is, of the order of the curvature scale, the so-called Hubble radius. The 
curvature scale is the characteristic length and will be also used as unit of length. For a scale 
factor dependence aft) oc t n in the time interval (t,t max ) the event horizon reads for t -C t max 


devent (t, tn 




H~\t) (h^) 1 n , 


for n < 1 
for n > 1 


( 1 . 2 ) 


where H = d/a = n/t is the Hubble scale. For an accelerated expansion, n > 1, the event 
horizon is roughly d even t(t) ~ H~ l (t). Obviously, the dependence is also similar for the a(t) oc 


- 3 - 











Figure 2: The figure illustrates schematically the initial conditions problem for the plateau infla¬ 
tionary models as the Starobinsky R 2 where F[nf ~ 10 ~ w Mp. The delayed inflationary dynamics 
imply that right after the Planck time hundreds of billions of causally disconnected regions (CDR) 
have to be homogeneous in order not to spoil the onset of inflation (Dh = -Dhomog(tp))- 


e Ht case. The importance of the event horizon is that it protects the initial smooth patch from 
the outside inhomogeneous regions where the gradients of the field are nonzero. Otherwise, 
if the event horizon had been unbounded, the inhomogeneities would have propagated and 
infested the initial smooth patch, rendering it unable to accommodate the inflationary phase. 

It has been actually shown that homogeneity on super-Hubble scales is required in order 
for inflation to start [15, 18, 19]. Indeed the inhomogeneity due to the gradients of the fields 
are critical and can prevent inflation. Inflation starts when V(< p) > ^(V(^/a) 2 ~ ( 6ip/aL ) 2 
where L the comoving wavelength of the homogeneity and 5p the typical change in ip. At the 
onset of inflation it is y/3HMp ~ y/V > 6ip/aL hence 


aL 

W 1 


> 


5ip 

3Mp 


(1.3) 


This qualitative estimation demonstrates that field variations 5<p > Mp cannot have a wave¬ 
length smaller than the Hubble radius. Any such inhomogeneity has to have a wavelength 
A p h = aL > few x H = £ H -1 . For inflation to occur the homogeneity has to be assumed 
on super-Hubble scales [15, 18, 19]. 

If inflation is unable to start at energies close to the Planck scale, as it happens at the 
Starobinsky R 2 model, then the minimum size of the initial homogeneous patch has to be 
much larger than Ip. The Starobinsky potential has two regions, a bounded from above 
plateau region (ip ^$> 1) that drives inflation at energy densities V(nf <C Mp, and a divergent 
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part (ip < 0) which, however, cannot constrain the kinetic energy of the held and no inflation 
takes place. Given the smallness of the Starobinsky inhationary plateau energy density, VJnF) 
one has to assume that a kinetic-energy domination regime preceded the inhationary phase. 
In such a case, V(ip) <C \ip 2 ~ ptotj the scale factor grows like t 1 / 3 until the domination of the 
plateau potential yielding an event horizon of size devent (iinit ~ ip, imax > Hnf) ~ 3 x 10 ^Hp 1 1 
where Hp 1 = \/3 Ip is the Hubble radius at the Planck time. Hence, one has to expel the 
density inhomogeneities at least 10 3 Hubble scales farther if the Universe has emerged from the 
Planck densities, Mp = (2.4 x 10 18 GeV) 4 . In particular the minimum required homogeneous 
region is 


-Dhomog(iinit) — d eV ent (iinit, Unf) “k [devent (ilNF ? ^max ) + H 1 (Unf)] 


a(t, 


init ) 


a(tiNF ) 


(1.4) 


or 


-^homog(iinit) — d eV ent(hnitUmax) ~k H (Unf) /, 


(1.5) 


where i max Unf a time deep inside the inhationary era. That is, inflation requires at Hnf 
a homogeneous patch of minimum radius 


£H 1 (Hnf) 


devent (j INF, f max ) + H (Hnf) 


( 1 . 6 ) 


which can exist only if the primary patch at ti n j t , £ H~ l (t inf) «(iinit)/o(UAfp)) is surrounded by 
a supplementary homogeneous shell of width equal to the event horizon distance d ev ent (^init > Hnf)> 
see Fig 3. The parameter £ is of the order 0(1). In particular in [15] the £ is evaluated to 
be £ > 3/2 for exponential inflation, i.e with equation of state w = —1. Here we find £ ~ 2 
for inflation to start at low energies, where Unf is the time that the equation of state of the 
helds drops below w = —1/3. 

For hat space and f; n j t ~ tp and tjNp ~ ^nf corresponding initially homogeneous 
volume, (4/3)7t H^ omog (tp), is at least 10 11 times bigger than (4/3)7r/p 3 which means that, 
initially, hundreds of billions of causally disconnected regions were much similar without any 
dynamical reason. Briefly we call them Causally Disconnected Regions (CDR). We consider 
the Ip as the causal horizon at Planck times. For an open Universe the number of CDR 
required to be homogeneous is even larger while for a closed Universe the number is decreased 
about an order of magnitude, albeit the CDR remains formidable large. In fact, these are 
much special initial conditions for the R 2 model and any similar plateau potential inhationary 
models. 

On the same footing with the CMB homogeneity reasoning, one concludes, for a trivial 
topology, that it is respectively unlikely that a Universe began in a chaotic state and has 
reached the homogeneous state required for the plateau inflation to start. 


1.2 Outline of the objectives and results 

One of the radical implications of the CMB data interpretation is that the R + R 2 success¬ 
fully ht PLANCK results [6] hence, it may provide insight into the effective description of the 
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homogeneous 
initial patch 



Figure 3: A schematic illustration of the minimum homogeneous patch at two different times: 
right after the Planck time, t = tp, and at the onset of inflation, £inf ~ 10 5 fp. Its radius at tp is 
Dh = P^homog {tp) and at fiNF is £-^inf ^ 


fundamental theory of gravity at the particular energy scales. In this paper we focus on this 
possibility examining the R + R 2 pure gravity and pure supergravity [20-22] , and their status 
related to the initial conditions problem. Indeed, supergravity, as the low energy limit of string 
theory, is a well-motivated effective description of gravitation for energies relevant to inflation. 
Our calculations are performed in the Einstein frame. 

It would be rather straightforward to invoke couplings to matter and build by hand scalar 
potentials which exhibit initially a chaotic inflationary phase (or inflation of the old type) and 
subsequently the final plateau inflationary phase; or to invoke non-minimal couplings that 
yield an effective potential which is flat enough to fit the data with inflation starting from 
Planck densities [23-25]. However, by turning to such models one abandons the simplicity and 
universality of pure gravitational (or supergravitational) microscopic description of inflation. 
In this work, we focus exclusively on the pure R + R 2 (super)gravity models. 

Our goal is to revisit the problem of the initial conditions in a supergravity setup, and 
in the case that the spatial geometry of the primal universe has not been flat. Despite the 
observed flatness of the present Universe we should expect that the initial patch had a non-flat 
geometry before the onset of inflation. It is inflation itself that justifies the post-inflationary 
flatness. We consider open and closed FLRW spatial background geometries. Other sugges¬ 
tions for non-trivial pre-inflationary topologies can also be found in the literature, see e.g. 
[26, 27] for compact flat or open Universe. 

Minimal supergravity has two different formulations: the old-minimal [28, 29] and the new- 
minimal [30]. One of the objectives is to examine whether the embedding of the Starobinsky 
model in minimal supergravity renders it more motivated in terms of the assumptions usually 
requested for the inflationary initial conditions. We report an affirmative answer to this 
question: the initial conditions are significantly relaxed, however, not fully addressed. The 
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reason is the presence of the dynamical pure supergravitational “auxiliary” 3 fields. 

Let us explain the origin of these fields. In supersymmetric theories there exists a class 
of bosonic or fermionic fields which do not propagate and are in principle integrated out. 
Their existence is essential for the off-shell closure of the algebra and for the construction 
of generic couplings. In fact superspace methods by construction give rise to these fields. 
Supergravity theories, being supersymmetric, also include auxiliary fields, and the difference 
in the auxiliary held sector is the root of the difference between the old-minimal and the 
new-minimal formulations. When higher curvature terms are introduced, the auxiliary fields 
might pick up kinematic terms. In fact, for R + R 2 supergravity this is exactly what happens. 
Therefore, when building R + R 2 theories of old-minimal or new-minimal supergravity, by 
construction one will find additional propagating degrees of freedom; these are a by-product of 
supersymmetry and impossible to avoid. In this work, instead of treating these fields as merely 
an eccentricity of supersymmetry, we show that they have important physical consequences. 
Moreover, since we are working with supersymmetric theories, these fields have to reside inside 
appropriate supermultiplets. Indeed, these supermultiplets can be uncovered by turning to the 
dual description of R + R 2 supergravity, which is standard supergravity coupled to additional 
matter fields: 2 chiral multiplets for the case of old-minimal supergravity [31-36] or one massive 
vector multiplet for the case of the new-minimal [34, 37, 38]. To perform the analysis in the 
following sections we will employ the dual description of the R + R 2 supergravity theories, 
but one should bear in mind that these additional superfields have a pure supergravity origin 
and therefore no additional sector is invoked other than pure supergravity. 

Turning to the dynamics of these fields, although they cannot initially drive inflation, 
they do implement a relatively fast expansion rate. Inflation starts naturally after a period 
of “kinetic-potential energy balance” generically yielding much more than 60 e-foldings. The 
relaxation of the initial conditions is greater in the old-minimal embedding where the dynamics 
of the “auxiliary” fields is described by scalar bosons. In the new-minimal, the extra bosonic 
fields include a gauge held that leads to an anisotropic expansion which, though it does not 
prevent inflation from starting, it ameliorates the initial conditions problem to a smaller extent. 
We also note that although the Universe can start with a potential energy V{ip) not too much 
smaller than the Planck energy density there may be no eternal inflation in the supergravity 
Starobinsky model. 

A homogeneous FRLW patch features either an open or a closed geometry. For the later 
case, the conditions for initial homogeneity are translated to conditions for the initial size of 
the post-Planckian closed Universe because the curvature term has to be always subdominant 
in order the collapse to be avoided. For the former case, however, the curvature term can 
dominate. In such a case the radius of initial homogeneous region though remarkably small 
the number of the CDR is found to be large again. 


3 Actually, the adjective “auxiliary” is literally wrong. In R + R 1 supergravity these fields propagate and are 
not auxiliary at all. However, we keep this term using quotation in order to keep the standard supersymmetry 
terminology. 
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Figure 4: The shape of the scalar potential for the old-minimal R + R 2 supergravity Starobinsky 
model. The one direction in field space corresponds to the inflaton <p and the second to the 6-field 
that makes Planck densities accessible. The potential is depicted in Planck units, Mp = 1. 

In the next section we present the pre-inflationary dynamics of the R + R 2 theory embedded 
in the old-minimal supergravity framework and specify the initial conditions for inflation to 
start. In section 3 we repeat the analysis for the new-minimal supergavity case where the 
pre-inflation expansion is anisotropic. Closed and open background FLRW geometries are 
considered in section 4 and last, in section 5, we conclude. 

2 Old-minimal R + R 2 supergravity: the F sugra ij 2 

The Planck data [6] favour inflationary models which predict a small amount of gravitational 
waves (small r) and spectral index n s = 0.96. This has given rise to the concern [9] that 
inflation is not an invulnerable candidate to solve the fine tuning problems of the Big Bang 
since it may have a fine tuning problem of its own. This is related to the fact that models which 
predict small r, for example the Starobinsky model of inflation [17] (see [39] for different models 
with similar effective description during inflation) have a flat potential Vp 2 that predicts a 
tensor-to-scalar ratio r* = 0.0033 at the pivot scale. For this value of r, the energy density 
of the plateau can be found from the relation F* = 1.5ir 2 A s r*Mp ~ (0.8 x 10 16 ) 4 GeV 4 
~ 1.2 x 10~ 10 A7p. This gives rise to a cut-off much smaller than the Planck mass 

Finf < 1.2 x 10 -10 Mp. (2.1) 

The natural way to overcome this is by somehow being able to have F ~ Mp at the initial stage 
of the Universe. It is easy to see that in the Starobinsky model (1.1) Planck scale potential 
energy densities can never be realized. This leads to a required fine tuning of the inflationary 
initial conditions. A different approach has been suggested in [40]. In the following we will 
show that supergravity offers a natural relaxation to this problem. 




The old-minimal supergravity multiplet [28, 29] contains the graviton (e^), the gravitino 
{ipm), and a pair of auxiliary fields: the complex scalar M and the real vector b m . As we 
have explained, supersymmetric Lagrangians with curvature higher derivatives also introduce 
kinematic terms for the “auxiliary” fields M and b m . The embedding of the Starobinsky 
model of inflation in old-minimal supergravity in a superspace (see for example [20]) approach 
consists of reproducing the Lagrangian 


Mi 


i vi M ' 2 


( 2 . 2 ) 


This is achieved by [32-36] 


C = -3 Ml 


d A 0 E 


i - A,nn + ^iz 2 iz 2 




3m 4 


(2.3) 


Modifications and further properties can be found in [41-51]. Note that R + R 2 non-minimal 
(20/20) supergravity has been constructed in the linearized level [53], using a formalism which 
can also describe higher superspin theories [52, 54, 55]. Lagrangian (2.3) when expanded to 
components yields R 2 terms and kinematic terms for M and b m . One may work directly 
with (2.3) but it is more convenient to turn to the dual description in terms of two chiral 
superfields: T and S. In other words, the superspace Lagrangian (2.3) has a classically 
equivalent description as standard supergravity coupled to additional superfields [32], The 
equivalent description of the above higher curvature supergravity reads 


£ 


d 2 0 2 £ 




+ c.c. + 


d 2 @ 2£W + c.c., 


(2.4) 


with Kahler potential 

K. 


and superpotential 


—3 Mp In 


T + T SS 1 5 2 5 2 ] 

+ M P M 2 + r M 4 J ’ 


(2.5) 


W = 6 mTS. 


( 2 . 6 ) 


During inflation the universe undergoes a quasi de Sitter phase which implies that supersym¬ 
metry is broken. In principle the identification the goldstino supermultiplet even though it 
is plausible is not always straightforward. An inspection of the properties of the model (2.4) 
during inflation shows that the goldstino multiplet is in fact the multiplet S. Moreover since 
the mass of the sgoldstino becomes large it can be integrated out [56, 57], leading to a non¬ 
linear realization of supersymmetry during inflation as was proposed in [58] . This amounts to 
setting [58] S = Xnl where = 0. Of course this effective description breaks down at the 
end of inflation. This idea has been further developed in [59, 60]. Note that the imaginary 
component of T is not integrated out due to the non-linear realization. It is strongly stabilized 
during the Starobinsky inflationary phase and therefore does not interfere with the dynamics. 
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Figure 5: The evolution of the supergravitational scalar fields ip (left panel) and 6 (right panel) for 
different initial values and F(c^i n it, 6; n it) = /Okin.init = 0.5ptot,init = Mp. The first field that evolves 
is the p while the b remains nearly frozen having an evolution timescale much larger. Once inflation 
begins around t ~ 10 5 fp the ypfield slowly rolls and the kinetic energy of the 6-field gets diluted. 
Initial values for the ip about zero can yield many more that 60-efoldings with ^>inf > 10 Mp and 
reach the minimum (<p,b) = (0,0) after about 10 lo tp. 


Eventually one finds the effective model 

e _1 C = — ^dipdip — ^m 2 Mp ^1 — e \Zi v/ Mp \ ^ ( 2 . 7 ) 

which is illustrated in Fig. 1. We remind that in order to find the Lagrangian (2.7) from 
(2.4) one has to set S = 0 and ImT = 0; for these values the fields are strongly stabilized 
[33, 34, 36]. From the Planck data [6] we get 

m ~ 1.3 x 10 -5 Mp. (2.8) 

Now we want the potential energy to become V ~ Mp for appropriate field values. The 
potential in (2.7) is unable to do this as we explained. 

In recent work [36] there has been found a new class of R-symmetry violating R + R 2 models 
which can both provide an inflationary sector and a hidden supersymmetry breaking sector, 
without invoking any matter superfields. The new properties of these models which distin¬ 
guish them from the R-symmetric R + R 2 old-minimal supergravity is that at the end of 
inflation the S field contribution starts to become important and the field configuration is 
driven towards the supersymmetry breaking vacuum. For these models it is also expected 
that the initial conditions problem is similar to the R-symmetric case that we analyse here. 

2.1 Evolution 

Now we allow for the ImT to take large values. First, for large ImT values this component 
also becomes dynamical and after the redefinition 

ip = 

b = V 6 ImT, 


- M P In [1 + 2 (ReT/Mp)] 


(2.9) 


10 











we have 


where 


e 1 C= -^f-R - \dipdy - \e 2 Vw /Mp dbdb - V sugraR 2 (tp, b ), 
^sugratfOP. b) = ^m 2 AfJ, ^1 - e ~^ v/Mp ^j + ^m 2 e~ 2 ^ /Mp b 2 . 


( 2 . 10 ) 


( 2 . 11 ) 


Second, the field S remains strongly stabilized and will not affect the evolution. Indeed for 
the mass of the complex held, S, we have 


m| = 


- 2,/ 2 

e V 3^771 

144 


2b 1 ( 24 + eV ) + 3 ( 24 + e VOip ( - 2e"V 3^(12 + Q + e V 5^(48 + Q 




( 2 . 12 ) 


which implies that for a moderately large value of the ((-parameter [33, 36] and for values for 
<p and b that give Planck-scale energy densities (2.21) the formula for the S-mass gives 

m 2 s > H 2 . (2.13) 


For the rest of our discussion we always assume that when the system evolves, the S held is 
heavy (2.12), stabilized at S = 0, and we do not take it into account for the system evolution. 
This is a rather typical setup in the old-minimal R + R 2 supergravity dynamics [33, 34]. 

We now turn to a hat FLRW background and study the evolution of the helds and of the 
spacetime. For 


2 5 (y/~g C ma t) 
V~9 S 9 mn 


the time-time component of the Einstein equation gives 


(2.14) 


3 H 2 M 2 p 


+ 7y e 2 ^ v/Mp b 2 + V sugraR 2 (99, b ) . 


Extremizing the action we take 


5C 


V=~9 


dn 


V^~9 


5C 


b{d m (t>) 


= 0, 


(2.15) 


(2.16) 


where cf> = ip,b and \J—g = a 3 . The equations of motion for the helds (p and b read 

<p + 3H<p + dVsngr aR2(y ’ 6) = — l e~ ^ <p/Mp b 2 , (2.17) 

dip V 3 

b+(3H - 2^| (pMp^b + e 2 ^ v/Mp h \ = q , ( 2.18) 
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and in a more analytic form 


and 


p + 3Hp + J-m 2 Mpe fl — e 

V 2 ' ' (2.19) 

_ y| e - 2 4^ (rn 2 b 2 - b 2 ) = 0, 


b + 3 Hb — 2 \j -pb + m b = 0. 


( 2 . 20 ) 


We want to study the evolution of this system which has been also discussed in [61-63] and 
in a different context in [64], 

We choose initial conditions such that the initial energy density, = Mp, is equally 
partitioned between the kinetic and the potential terms, 


^(^initj frinit) — /4dn, init — ^ Mp. 

An example of such a set of values that realize these initial conditions is, 


A/ 2 ■ 1 

(finite hnit) = (0, ) and firdt = hnil = / -1/ / ’ • 

m V 2 


( 2 . 21 ) 


( 2 . 22 ) 


The equations (2.15), (2.19) and (2.20) can be solved numerically. For the initial condi¬ 
tions (2.22) we have found exact numerical solutions, which are illustrated in Fig. 5. It is 
easy to see (also from the potential) that the f held, due to the large b values, has locally a 
run-away potential and starts to grow. Later, the b held starts to roll down its potential. This 
is the imaginary Starobinsky phase [61-63]. At some point the b becomes small and then the 
<p held stops increasing and we have the inhationary initial values 


(VhNF, &inf) ~ (13Mp,0) 


(2.23) 


for the (2.22) initial Planck-scale energy densities. Then a Starobinsky inhationary phase 
starts and is described by the Lagrangian (2.7), which naturally lasts for much more than 60 
e-foldings 4 . 

To summarize, the dynamical evolution the pre-inhationary stage consists of two phases: 

1. From Fjugrai? 2 — Afp to K. ugr£L R2 > m 2 Mp, both p and 6 participate in the evolution 
with the p held rolling hrst. 

2. At 14 U g r£L R2 ~ m 2 Adp starts the standard Starobinsky inhationary phase with p driving 
inflation, and b now strongly stabilized and integrated out. 

4 We note that there is no fine tuning problem here similar to that of hybrid inflation where special initial 
values for the fields are required in order N > 60 e-foldings to be achieved [65-68]. 
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Evolution of the equation of state 


Evolution of the scale factor 




Figure 6: The left panel shows the equation of state, w, for the supergravitational system of fields. 
The initial conditions chosen are equipartition of energy between kinetic and potential thus w = 0 
initially. For some period it is w ~ —1/9 and at fjNp — 0.7 x 10 5 fp the equation of state becomes 
w ~ —1/3 and acceleration starts. The nearly de-Sitter phase w ~ — 1 starts after 3 x 10 5 tp. The right 
panel shows the evolutions of the cosmological scale factor. The solid lines correspond to solutions for 
the scale factor of the conventional Starobinsky (lower, blue) and the Starobinsky supergravity (upper, 
red). The initiation of the accelerating phase is apparent after tiNF- The dashed lines close to the solid 
ones are the constant equation of state approximations. The lower blue dashed corresponds to the 
w = 1 equation of state and describes exactly the evolution of the scale factor before inflation for the 
Vr 2 case; the green dashed corresponds to constant w = —1/9 which approximates well the Esugra r 2 
case until, roughly, the onset of inflation. The upper black dashed line corresponds to background 
spatial geometry of negative curvature for the conventional Starobinsky plateau inflationary potential. 


2.1.1 Dynamics of the expansion 


The local conservation of the energy-momentum, X7 m T mn = 0, gives the evolution of the 
energy density which in an FLRW background reads p + 3 H{p + p) =0. When the energy 
density is dominated by a fluid with constant equation of state, w = p/p, the energy density 
changes with the cosmological scale factor as p(a) = Pinit( a / a mit) — From the Friedmann 
equation we take 


a(t) — Q-init 




-i 2/[3(l+m)] 


(2.24) 


where a(fmit) = «init, p(ainit) = Pinit for tj n i t < t. When the small constant term in the 
brackets is negligible, we see that the scale factor changes as a{t) oc and the energy 

density as 

4 Ml o , 

P - , P , 2 t 2 , for t > fjnit • (2.25) 

3(1 + w) z 

The time that signals the onset of inflation is set by the energy density of the plateau. 
When the total energy density is p ~ Vinf then inflation starts and from eq. (2.25) we take 
that 


flNF 


4 Ml 
3(1 + w) 2 


V T 


-1 

INF 


1/2 


lO^fp 1 = I0°ip. 


(2.26) 
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The precise time is actually smaller because acceleration starts when w < —1/3 that is when 
the kinetic energy density is half the potential energy density, hence pinf > Vinf and Hnf < 
10 5 fp. We find numerically that for the R + R 2 supergravity theory inflation starts when 

t mF ~ 0.74 x 10 5 M~\ (2.27) 


for equipartitioned initial energy densities. 

The system of the fields (<p, b) starts from nonzero values such that V (pinit,- frinit) ~ Mp. 
Due to the small mass of the b field, mb H , the cp will roll down the potential (2.11) which 
for constant b has the form 

V sugiaR 2&) ~ V 0 e~ 2 ^ /Mp . (2.28) 

Initially, the approximation 6 ~ 0 is a good one according to the numerical results. For the 
potential (2.28) the Friedmann equation (2.15) and the equation of motion (2.17) have an 
exact solution of power law form given by the expressions [69] 


a oc ai n itt n , n = 3/4 



(2.29) 


According to (2.24) we see that n = 2/[3(1 + w)} = 3/4. This corresponds to a barotropic 
fluid with equation of state w = —1/9, that is a negative pressure. Numerically we find that 
in the pre-inflation period the energy density decreases with a slower rate, 


w. 


sugrai? 2 1$ 1/9) ®sugrai? 2 (l) ^ ^ i Psugrai? 2 Pinit® ^ 


(2.30) 


because the actual system is a two-field one. 


2.1.2 Selfreproduction 

At this point we would like to briefly comment on the initial conditions for the R + R 2 
(super)gravity model that lead to the eternal process of self reproduction. After tjNF the <p- 
field rolls pretty slowly, as depicted in the equation of state w at Fig. 6, and the average 
quantum fluctuation |<5<p| ~ H/2ir ~ m/(47r) inside the domain of radius H can be larger 
than the classical variation of the inflaton |Atp| ~ ^3/2 Mp/N(ip) = 2^j2/3 Mp 
during a Hubble time, where N(<p) is the number of e-foldings before the end of inflation. 
Eternal inflation takes place for |<5cp| > |A<p| or equivalently for sufficiently large </?inf values 

<Pinf > 17.5 M P = <pi NF , (2.31) 

where the value (2.8) for m has been plugged in. For example, when the initial values of the 
field is </>init 0 and Vinit 1 then inflation never enters the self-reproduction regime. As it is 
listed in the table 1, in the old-minimal embedding, only values larger than <pj n it > 13.5 lead 
to the self-reproduction regardless the fact that initially the potential energy density is of the 
order of the Planck scale. 
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The event horizons 



Figure 7: The figures show the event horizon distances for the Starobinsky (blue) and supergravity 
Starobinsky (red) model. In the left panel, in a logarithmic plot, we see that the event horizon for the 
supergravity is about two orders of magnitude less than the conventional Starobinsky and it becomes 
nearly constant after t/jvF', this is evident from the dashed line that corresponds to w = 1 and describes 
the pre-inflationary evolution of the Starobinsky model. The right panel shows the event horizon for 
the supergravity Starobinsky for different initial values for the fields. When the p tot , init is entirely 
potential energy then the event horizon distance is smaller than in the equipartition case. The lower 
black dashed line shows the event horizon distance when the negative spatial curvature dominates. 
The event horizon distance is in Hp 1 = \/?>lp units. 


2.2 Initial conditions 


The standard lore of inflation is that it started right after the Planck era and magnified the 
volume of the Universe e 5N times, where N > 60 the number of e-foldings. 

If inflation begins at times t \nf tp, i.e. energy densities V <C Mp, as it happens with 
the plateau potentials, then in the period between the Planck time and t \nf the equation of 
state of matter in the Universe has to be w > —1/3. If w is nearly constant then the scale 
factor increases like a(t) oc t n where n = 2/(3 + 3 w) < 1. Then the part of the space that will 
be in causal contact until inflation begins has radius 


devent (^init i t INF i w) — Gr(tinit) 


UlNF 


n 


dt 


a[t ) 1 — n 


H-\t init) 


( Unf 
V ^init 


1 —n 


(2.32) 


for fjnit <S tjNF- 

2.2.1 R + R 2 

In the non-supersymmetric Starobinsky, the energy density has to be dominated by the kinetic 
term c/> 2 /2 as long as p > Vinf which translates into an equation of state w = 1. The same is 
true for all the similar single-field plateau models where U[nf r\j 10 10 Mp. In such cases the 
energy density falls fast as p oc a -6 while the expansion is especially small a(t) oc t 1 / 3 . This 
implies that the initial patch is sensitive to inhomogeneities far away because the intervening 
space expands too slowly. In particular the event horizon radius, when we integrate from the 
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Planck time till the onset of inflation, Hnf — 0.6 x 10 5 t P , is 

1 ft \ 2 / 3 

devent{tmit = tp, HnF, ^ 2 ) ^ ~H~ l {t P ) ( J ~ 766 H p l (2.33) 

where Hp 1 = \/3 l P . Actually, the event horizon continues increasing after the initiation of 
the accelerating phase, tiNF, though much more slowly, and does not get a constant value 
until the equation of state becomes w = — 1, see Fig. 6 and 7. We find numerically that the 
total event horizon radius at the Planck time is 


devent(^init — t Pl t max , WR 2 ) — 2820 Hp — 4884Zp , 


(2.34) 


where t max Hnf a moment inside the inflationary era. Hence, the minimum initially 
homogeneous region required for inflation to start has radius 

-^homog(^P) ^ R 2 ) = d even t(t P , tmax) 4“ H (Hnf)^ 

~ 4884 Zp + H 1 (Hnf) ~ 8.7 

where g/Rnf) ~ 47 a{t P ) and ~ 3 Hnf- The minimum number of the causally 

disconnected regions (CDR) required to be homogeneous is 


(tp) 

RNF ) 

x 10 3 1 p 


(2.35) 


1(1 at ( T^homog; p2 ) 

Cfl at (Zp) 


-7tD 3 

3 4 ho 3 mog ~ (8.7 x 10 3 ) 3 ~ 7 x 10 11 CDR, 

gP/p 


(2.36) 


which manifests the initial condition problem for the plateau potentials such as the Starobinsky 
R 2 model. 


2.2.2 R + R 2 supergravity 

On the other hand, in the Starobinsky supergravity model the pre-inflation expansion of 
space is much faster (2.30) and the event horizon is much smaller. An approximation (2.29) 
is to consider a constant equation of state w = —1/9 which yields an event horizon radius 
devent oc (fmax/ip) 1 ^ 4 ~ 87 l P for t max = 10 5 fp. An exact result can be obtained numerically 
for the varying equation of state rc sugra p 2 (2.30). When we integrate from the Planck time 
until the beginning of inflation, which is found to be Rnf = 0.74 x 10 5 t P , we numerically take 

^event (^init = tp, h nf, W su graR2 )~29 Hp 1 . (2.37) 

As in the R + R 2 gravity case, the event horizon increases as long as w > — 1. It remains 
constant when the field configuration lies in the plateau with vanishing kinetic energy, w = — 1, 
see Fig. 6 and 7. The numerical value of the total event horizon reads 

devent (bnit — tp , f maX! ^sugraP 2 ) — 46 Hp — 80 l P . (2.38) 
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The minimum initially homogeneous region required for inflation to start in the supergravity 
case has radius 


Dhomog(tPj W sugrai? 2 ) — devent (tp , fmax )+H 1 (Rnf) u P \ 

a(t INF ) 

~ 80 Ip + H 1 (Rnf) 


(2.39) 


5407 


98 l P 


where a(fiNF) 5407a(tp) and H ^inf) 4/3tiNF- That is, right after the Planck time 
the initial homogeneous volume is required to have radius at least 68 times the Planck length. 
The minimum number of the CDR is here 


!^at(T?homog) ^sugraR 2 ) _ S^^homog 


y&SLt(lp) 


|vr/3 


(98) 3 ~ 10 6 CDR. 


(2.40) 


Compared to the non-supersymmetric case, in the R + R 2 supergravity the required initial 
homogeneous volume is about half a million times smaller, 


#CDRp2 _ -Pg omog (f P , w R 2) 
#CDR sugra p 2 d)^ omog (tp, w 

sugraR 2 ) 


(2.41) 


To outline, in the R + R 2 supergravity the initial conditions problem though significantly 
ameliorated (about one million times) it persists. The evolution of the event horizon for the 
Starobinsky supergravity as a function of time and the initial conditions (the energy partition 
between kinetic and potential) can be seen in Fig. 7. 


3 New-minimal R + R 2 supergravity: the VnsugraR 2 


The new-minimal supergravity multiplet [30] contains the graviton held e^, the gravitino 
which are physical fields, a real auxiliary vector A m which gauges the 17(1) R-symmetry 
and a two-form auxiliary held B mn . The two-form appears here only through the dual of 
its held strength H m . In the theory with no curvature higher derivatives A m and H m are 
integrated out. This does not happen when we introduce higher curvature terms. In this case 
a combination of these helds becomes propagating [37]. The Starobinsky model of inflation in 
new-minimal supergravity reads in superspace [34, 37, 38] 


2 M 2 p J d 4 6EV R + | J d 2 d£W 2 (V R ) + c.c., 


C = - 


and in component form for the bosonic sector we find 


(3.1) 




1 


C = M 2 p (-^R + 2 A a H a - 3 H a H a )+^(-R + 6 H 2 ) 1 - jF 2 (A~), 


r 2\2 


a 

4 


(3.2) 


for A m = A rn — 3 H m . It is easy to verify from (3.2) that when there are no curvature 
higher derivatives present, e.g. the limit a —> 0, the helds A m and H rn vanish on-shell. The 
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(a) 

Figure 8: The evolution of the supergravitational fields ip (left panel) and A z (right panel) in the 
new-minimal embedding for different initial values and F(<^i n it,-Az,imt) = Pkin = 0.5ptot,init = Mp. 
The vector field A z gets redshifted from the expansion and the scalar <p rolls a smaller distance in 
field space than in the old-minimal case ruling out initial values ^i n ;t < 1.5 Mp because not enough 
e-foldings are achieved. Once inflation begins around t ~ 10 5 tp the (/ 5 -field slowly rolls and the vector 
field A z gets completely diluted. 



Lagrangian (3.2) is classically equivalent to a Lagrangian where no higher derivatives are 
present, in particular for a = gp- we find 

e~ l C = -^R - \f\V) - \dipdip - ^ (l - e-^ 2 - 3 g 2 e~ 2 ^ V m V m , (3.3) 

where V m = A m — 3 H m . This is standard supergravity coupled to a massive vector multiplet. 
Various modifications can be found in [70-73]. We will use the dual description (3.3) for our 
discussion of the R + R 2 supergravity. 

3.1 Evolution 

From the Lagrangian density (3.3) we see that the only way to increase the energy density to 
Mp is by giving an initially large value to the vector V m . In this scenario we choose the gauge 

Vo = 0, (3.4) 

and we take the ^-spatial axe parallel to the direction of the vector 


K = A z (t)6f. (3.5) 

By giving to the vector a non-vanishing value a direction is singled out from the other two 
perpendicular in the spatial space. This implies that the metric will be described by two scale 
factors 


ds 2 = —dt 2 + a 2 (t)[dx 2 + dy 2 ] + c 2 (t)dz 2 , (3.6) 

hence, an anisotropy is created. Here we have identified this direction with the z-axis. 
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The Einstein equations read 


with 


a\ ac 1 
a) +2 “" = 


a c Mp 1 


c a ac 1 1 

- + - + — = -— Px = -— Py , 


c a ac 
2 


. a 


a 


2 —hi- 


1 

M 2 


Pz, 


P= 2^ + 2 


1 / Az 

c 


T ^lsugrai? 2 (Ti Az ) j 


1.2 1 (A 


Px Py cJp T 14isugrai ? 2 {jpi Az)-, 


(3.7) 

(3.8) 

(3.9) 

(3.10) 

(3.11) 


Pz = 7i<P 2 -7i I — I - *9* ( 1 - e 


1 / A 


+ 3, 2 e- 2 4^(^) . (3.12) 


The energy densities and pressures are defined by the energy-momentum tensor Too = P- 
T xx = Px/a 2 , Tyy = py/a 2 , T zz = p z /c 2 for 


1 


1 


Tmn — 2 ^ mn ^ ^k^dlipJ + yF m iF n ^QmnF'klF 

— 9mn ^nsugrai? 2 {Pi •Az) + ^^nsugrai? 2 ^™^’ 


where 


Ksugra R 2 (PiA z ) — -g 1 - e 


+ 3 g e 


2„—2 


A z 

c 


(3.13) 


(3.14) 


and the b )( sugra ^2 (</h-4z) denotes the derivative with respect to the Lorentz invariant quantity 
V At V M = (A z /c) 2 . Extremizing the action, the held equations read 


We will set 


.. I r,® c\ ^UnsugraR 2 i}Pi A z ) 
v+ ^a + c) v + - 85 - = 0 - 

A z + | 2- - -) A> + 0V ™^«A' Az) = 0. 
a c) oA z 


g 2 = - m 2 ~ - (1.3 x 10 5 ) 2 Mj>. 
6 6 


(3.15) 

(3.16) 

(3.17) 


Starting from V(^imt, A z { n it) ~ Mp the hrst held to roll down the slope of the potential 
is the scalar ip. Initially the vector held has a small mass, m^ z <C H, and will stay nearly 
frozen. In this period the effective potential is 


Wa^M~Voe 2 ^ /Mp - 


(3.18) 
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Figure 9: The evolution of the scale factor (left panel) and the event horizon (right panel). The 
anisotropic expansion of the Universe is manifest which implies that the event horizon distance is 
respectively anisotropic. After the onset of inflation f INF the scale factors evolve similarly and the 
anisotropy gets diluted. The event horizon distance is in = y/3 Ip units. 


3.2 Initial conditions 


The evolution appears similar to the old-minimal case however, here, there is a background 
vector field with non-vanishing value which breaks the isotropy of the space and, as described 
in eq. (3.7) -(3.12), the scale factor in the directions parallel and perpendicular to the vector 
evolves differently. The expansion rate of each direction depends on the initial conditions. 
When pi n it — Unit 3> Pkin, init the z-component of the pressure, p z , is positive while the p x 
and p y are negative implying that the c(t) scale factor grows faster than the a(t) giving at the 
Universe a cigar-like shape. On the contrary, if pj n j t — pki n , init ^ Unit the ^-component of 
the pressure is negative while the p x and p y are positive and the a(t) scale factor grows faster 
than the c(t) and the shape of the Universe is a pancake-like one. We numerically solve the 
system of the equations and the evolution of the two scale factors, for pkin, init = (1/2) Unit 
can be seen in the Fig. 9. Accordingly, the event horizon distances change in the z-direction 
and the x — y plane. The volume of the event horizon is albeit found not to be much sensitive 
to the partition of the initial energy density. 

We can also redefine the two scale factors as a = e a+l3 and c = e a ~ 213 where the e“ can be 
seen as the isotropic scale factor and (3 the deviation from isotropy. The zero-zero component 
of the Einstein equation (3.7) can be recast into [74] 


3 Ml 3 


3Mp 


C P + Mn) 


(3.19) 


The 3 6 = (2a/a + c/c) = 3d represents the volume expansion rate, Vol/Vol, and o 2 = 3/3 2 
is the one half of the sum of the shear components, ( J t squared, where c,; = dj/ at — 6 . If the 
energy density of the vector field that feeds the anisotropy is vanishing then the background 
Einstein equations imply that /? oc (e“)“ 3 , hence the effective energy density of the anisotropy, 
PAN = Mpcr 2 , redshifts rather fast, similar to the kinetic energy density redshift. 
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The scalar field equation of motion reads 


ip + 36(f) + 


•^^nsugrail 2 ijPi *^z) 

dtp 


= 0 . 


(3.20) 


As we can see from the analogue of the Friedman equation (3.19) the Hubble friction has been 
replaced by the "volume expansion friction" which is larger than the usual value (p/SMp) 1 ^ 2 . 
The anisotropy aids the slow roll of the (p field, which at first sight may look somewhat helpful 
for the amelioration of the initial conditions problem. Although it would be rather interesting 
to find a slow roll regime due to the presence of the pan, the anisotropy does not yield any 
helpful increase in the expansion rate of the Universe. The numerical results are shown in 
Fig. 8 and 9. We finally mention that the primordial anisotropy described by the metric 
(3.6) corresponds to a Bianchi spacetime which does not prevent inflation from starting [15], 
a result which is also manifest in Fig 9. 

The minimum homogeneous region required at f; n it = tp for inflation to start at Unf in 
the new-minimal supergravity has volume 


(3.21) 


where devent the event horizon in the x — y plane and dgvent in the ^-direction. At the onset 
of inflation the contribution of the vector fields in the energy density is subdominant and the 
initial patch that gets inflated has a spherical volume V = |7riW 3 (iiNF). However, the initial 
homogeneous volume (3.21) is an oblate spheroid. 

The number of the causally disconnected regions is found numerically to be of the order 


14omog (t P ) ~ iq7 cdr 

P d P 


4 The curvature term 

An homogeneous and isotropic spacetime is described by the FLRW metric 
ds 2 = — dt 2 + a 2 (t ) + r 2 {d0 2 + sir^Odcj) 2 )^ , 

and the evolution of the scale factor is given by the Friedmann equation 

/a\ 2 p k 
\a/ 3 Mp a 2 


(4.1) 


(4.2) 
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Initial Values (old) 

^sugrai? 2 (^init? ^init) ^ 1 

Initial Values (new) 

^nsugrai? 2 (^init ?-^init) ~ 1 

N> 60 

Selfrepr. 

<^INF 

<^init > 4 

binit > 2 x 10 6 

ifinit > 13.5 
dinit > 5 X 10 9 

Yes 

Yes 

<PINF > 17-5 

- 8.5 < ipinit < 4 

10 2 < binit < 2 X 10 6 

1-5 < <Pinit < 13.5 

3 x 10 5 < d init < 5 X 10° 

Yes 

No 

5.5 < ipiNF < 17.5 

<Anit ^ —8.5 
binit > 10 2 

ifinit <1.5 
di„it > 3 X 10 5 

No 

No 

VINF < 5.5 


Table 1: The table presents the constraints on the (old- and new-minimal) supergravity scalar field 
space, in units Mp = 1, for sufficient number of e-foldings starting from Planck-scale initial energy 
densities for the potential and the kinetic terms. The field values that lead to a selfreproduction regime 
are also shown. 


When k ^ 0 the space is curved and the volume enclosed in sphere of radius d = a\0 differs 
from that of the flat space. It is [75] 



/ 

^closed 

= 7ra 3 (2xo - sin2xo), 

0 < Xo < 7T 


V = < 

^flat 

4 S 
|7ra , 

Xo = d com 

(4.3) 


^open = 

= 7ra 3 (sinh2xo - 2xo), 

0 < X o < °o • 



In the closed and the open case the xo is a dimensionless anglular coordinate and the scale fac¬ 
tor a is the dimensionful radial coordinate. In the flat space xo corresponds to a dimensionful 
comoving distance. For a definite xo it is 

Fclosed(,\o) < t4at(,\o) < KpenC\o) • (4-4) 

When xo 1 the volumes do not differ much while for xo > 1 the difference becomes 
important. 

If the energy density of the Universe is larger than the critical density, p CI n = 3 H 2 Mp, 
the Universe has a closed spatial geometry (k = 1), whereas for p < p cr i t an open ( k = — 1) 
one. In the special case that p = p cr ^ we have a flat geometry ( k = 0). In the preceding 
sections we have implicitly assumed that in the pre-inflationary stage it is | p — p C rit| 1 - 
This accounts for a rather special initial condition. Nevertheless the results obtained can be 
accordingly translated into the case that the spatial curvature is nonzero. We mention that 
the present data find no evidence for any departure from a spatially flat geometry [7] 

I^aI < 0.005, (4.5) 

where Q a = 1 — Ptot/Pciit- It is actually inflation itself that addresses the puzzle of the 
observed flatness of the Universe. The Ua decreases exponentially with time during the 
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accelerated expansion. In the models studied here it is very natural to expect many more 
than 60 e-foldings and this fact diminishes the chances for any observable deviation from 
flatness. However, before inflation a homogeneous initial patch is expected to feature either a 
closed or an open FLRW geometry. 


4.1 Closed Universe 


When the Universe has a positively curved geometry, the Friedmann equation can be written 
as 

3 MpH 2 = p - p c iosed, (4.6) 

where /Closed = 3 Alp/a 2 . At the moment fturn the Universe reaches its maximum size and 
a(tturn) = 0. There, the energy density has the value p(ttum) = 3AfJ,/at urn and the evolution 
turns from expansion to collapse. Inflation has to start before ft U rn, that is p(Unf) > Pclosed- 
Assuming a scaling of the energy density of the form 


P = Pinit 


Ginit 

a 


n 12 / t \~ 2 

~ n* \Mp) ’ 


(4.7) 


where the second equality follows from eq. (2.24) for t f; n it, we find the expression for the 
initial radius in terms of the aturn radius 


Ginit = ^3Mj,alJ n /n . (4.8) 

The turnover is postponed to the far future if p(Hnf) <1 3/aj( urn , i.e. a turn > (n/2)tiNF> that 
gives the relation between the time of inflation and the minimum initial radius, a(t = t; n ;t) = 
ai n it, of the closed Universe 

Ginit,min ~ (|) . (4.9) 

The above expression (4.9) is an approximate one because a subsequent accelerating phase 
takes over. The equation of state changes continuously from the value w = n/ 3 — 1 to the 
value w < —1/3 and the expression (4.8) is not exact when the scaling of the energy density 
deviates from (4.7). The (4.9) yields an overestimated value. Below we present the cases with 
and without supergravity separately. 


4.1.1 R + R 2 

When the inflationary energy density is bounded from above, as it is in the Starobinsky model, 
Uinf ^ 1.2 x 10 -10 Mp, then apparently, for our Universe to survive, the turnover energy 
density has to be in lower values: pturn < 4 ^nf — 1-2 x 10 _10 Mp. An approximate estimation 

is a turn > 0(^3Mp/Vinf)- The numerical estimation yields a turn > g(Hnf) — 1-3 x 10 5 Zp 
and Unf — 0.6 X 10 5 tp. Before inflation, higher energy densities are in the form of kinetic 
energy of the scalaron held, hence 
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For pkin,init = Mp and from (4.8) we take an estimation for the value of the initial radius of 

2/3 1/3 

the closed Universe aimt ~ a tum Ip ■ Numerically we find 

Omit > Qinit,,min ~ 2.5 X 10 3 Ip . (4-11) 


This is a rather large radius. It means that when the initial Universe emerged from the 
Planck, possibly quantum gravity, era it must have had a radius a(tp ) = oi n it of at least a few 
thousand times the fundamental Planck length. This is a formidable radius for theories that 
attempt a description of the quantum genesis of our Universe [2, 76, 77]. If it had any smaller 
size it would have collapsed before inflation begins. 

The minimal size of the initial radius a; n it depends on the time inflation starts. For a 
Universe dominated by the kinetic energy, the energy density falls like p oc aU 6 and from (4.9) 
we take for Hnf bnit = tp 


2/3 ,1/3 


®init,min(^init — tp,tj Nf) ~ ^INF ^ 


(4.12) 


We have implicitly assumed that the homogeneous volume is the entire volume of the 
closed Universe, 2-7T 2 af nit , because we considered an FLRW evolution for the whole spacetime 
not only for an initial patch. Here as well, the closed Universe with radius (4.11) contains 
billions of initially causally disconnected regions. The minimum volume (4.3) that can be 
considered to be in causal contact is the one enclosed in a sphere of radius d = aimt.minXo = Ip- 
It is xo = lp/ ®init,min ~ 1/2500 <C 1 hence 


Fclosed(®init X0 — Ip) — ^ Closed (oinit 71) — 27T^a,; n j t . 


(4.13) 


The entire volume V^ sed = Klosed(ainit,min tt) contains 


T/tot 

R closed _ imt,min -7 s , mlO rT»r> 

wsw~ 7 x 10 CDR 


(4.14) 


homogeneous causally disconnected regions (CDR) for /; n j t = tp as depicted in Fig. 10. 


4.1.2 R + R 2 supergravity 

In the Starobinsky supergravity model, potential energy density values I^ugrai? 2 (<p) ^ f^NF 
are possible 5 . We assume that the potential energy after the Planck era is of the order of 
the Planck mass to the fourth power, p- in a = Alp. An approximation for the scaling of the 
supergravity energy density is given by the expressions (2.29) that is 


Psugra ~ Pw= — 1 /9 — Pinit 


Qinit \ 8 / 3 
a ) 


(4.15) 


’In the supergravity case also holds Unf M %, however KugraflUv 5 ) ^ Uinf values do not rule out 
inflation as it happens in the conventional Starobinsky model. 
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Hence, ai n it ~ ( hmn ^p l - Asking again for p tum < Uinf — 1.2x 10~ 10 Mp we find an estimation 
for the minimum initial radius. The dependence of the aini^min on the time inflation starts 
has the following form for Hnf iinit = tp 

Oinit,min(^init — ^PUlNF) ~ ^INF ^ — 1/9. (4.16) 

The above accounts for a conservative approximation. Numerical estimations for the super- 
gravitational system yield aturn > o(Hnf) — 10 5 Ip, Hnf — 0.74 x 10' 5 tp and the lower bound 
on the minimum radius of the 3-sphere is 


®init ^ ®init,min — 1*7 X 10 Ip , for ^sugraP 2 ■ 


(4.17) 


We mention that when the spatial geometry is curved the scale factor corresponds to the 
radial coordinate and it is dimensionful. The ratio a(f)/a; n ;t denotes the number of times the 
radius of the closed Universe has increased compared to the initial radius. 

Compared to (4.11) the bound (4.17) on a; n it = a{tp ) accounts for more than one hundred 
times less severe condition for a closed Universe to reach the inflationary period before it starts 
collapsing. Again here, we have assumed that the entire space is an homogeneous three- 
dimensional sphere with thousands of initially disconnected regions. The minimum volume 
(4.3) that can be in causal contact is enclosed in a sphere of radius ttmit.minXo = Ip- It is 
Xo = Ip/ Omit,min ~ 1/17 and also here 

bclosed(Oimt X0 = Ip') — ^P^p ^closed (®init 7r) = 27T^a,; n j t . (4.18) 

The entire volume V^ sed = U c i ose d(oinit,min tt) contains 


T/"tot 

^closed 

^closed (^p) 


27r 2 a 3 


init,min 


|vr/3 


~ 2 x 10 4 CDR, 


(4.19) 


for fjnit = tp, as depicted in Fig. 10, that is few million times less CDR than the R 2 case. 


4.2 Open Universe 

The initial patch that has been inflated may locally resemble a geometry of negative curvature. 
In this case the Friedmann equation reads 

3 MpH~ = p + Popen j (4.20) 

where p 0 pen = 3 Mp/a 2 . Here, the space corresponds to a hyperbolic plane that has an infinite 
volume, though we are interested only in the local geometry of space inside the Hubble radius 
not globally. The scale factor can take small values without any fear that the Universe will 
collapse. The p open can dominate over the energy density of the Universe and dilute the matter 
till the energy is redshifted to the value of the inflationary plateau U[nf- Then the inflationary 
evolution takes over, spacetime becomes approximately de Sitter and the negative curvature 
term, represented by the p open , asymptotically vanishes. 


- 25 - 




This case yields a rather small event horizon distance. In particular the curvature term 
scales like 1/a 2 and the a(t) grows like t. In the eq. (1.2) the power is n = 1 and the event 
horizon untill inflation grows only logarithmically with time 


devent (Cnit , b n ax ) — o(bnit) 


^ ^init 


dt 

ait) 


— a(tinit) 


In 


( UNF 

V bnit 


+ 


' ilNF 


dt 

ct(t) _ 


(4.21) 


For i max > 10 5 tp tlNF we find numerically 


^event (bnit — tp, fmax > Hnf) ~ 12.5 a(fp) , 


(4.22) 


where a(b n it = tp) = a(tp) is the real value of the hyperbolic space radius. The minimum 
required homogeneous regions has radius 


F^homog(tp) — devent (tp, tmax ) + H 1 (tlNF) 


q(tp) 

a(tlNF ) 


14 a(tp ), 


(4.23) 


since a(t) ~ t. The (4.23) seems to be a comfortably small value. The imaginary radius of the 
open Universe can take small values such as the Planck length and expand fast and endlessly 
resulting, at first sight, in only a “logarithmic sensitivity" of the initial patch to the regions 
that can influence it. However, if the initial radius of curvature a(tp) is about the Planck 
length then the volume enclosed in a sphere of radius Hhomog ~ 14 a(tp) is remarkably large 
in such a highly curved hyperbolic space. Indeed it has to be a(tp)x o ~ 14 a(tp) hence the 
angular coordinate is \0 ~ 14 1 and the Euclidean space approximation breaks down. It is 


Uopen (a(tp)xo = 14a(tp)) = 7ra 3 (t P ) (sinh2\ 0 - 2xo)| X o=i4 — (7.2 x 10 11 ) ira 3 (t P ) (4.24) 


On the other hand, the minimum volume of a causally connected region is enclosed in a radius 
about Ip and has size 

Fopen (a(tp)xo ~ lp) = 7m 3 (tp)(sinh 2y 0 - 2x 0 )| i ~ 0.277ra 3 (tp). (4.25) 

\/3 

We took ai n it = a(fp) = lp\/ 3 in order that H 2 {tp) = l/(3/|>) and be in agreement with the 
flat and close Universe cases examined before where H 2 {tp) ~ Mp/3. Hence, we find 

V op en(D^ ) ~ ^ x 1Ql2 cdr ^ ( 4 . 26) 

Fopen (Jp) 

at ti n it = tp. This is a very large number. The (4.26) is actually larger than the CDR numbers 
the Starobinsky R 2 model yielded for the cases of flat (2.36) or closed (4.14) Universe. We note 
that in the closed Universe case the initial radius had to be large enough to evade collapse and, 
hence, the curvature term was subdominant. Here the space is hyperbolic and the Euclidean 
approximation cannot be applied for distances larger than the Planck length. 

If we consider a subdominant curvature term, p Qpen = 3 Mp/a 2 <C p then the flat space 
approximation is reliable and the results should be similar to those of flat and large radius 
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closed Universe. In particular, for the R + R 2 (super)gravity model, the -Dhomog should not be 
much different for k = 0, ±1. Then according to the hierarchy of the volumes in these three 
different geometries (4.4) we take again that 

# CDR(closed) < # CDR(flat) < # CDR(open). (4.27) 

Hence the hierarchy (4.27) holds either for a dominant or subdominant k = —1 curvature 
term. 

5 Conclusions 

In this work we studied the R + R 2 gravity and supergravity models for inflation (known also 
as Starobinsky models), which are characterized by a plateau inflationary potential and are 
particularly motivated after the release of the PLANCK 2013 results. However, they account 
for low energy scale inflaton models, requiring a rather extended acausal homogeneity in order 
for inflation to occur. We demonstrated that the problematic issue of the initial conditions is 
less severe if supergravity is realized in nature due to the extra directions in the field space that 
can implement a relatively fast expansion rate before inflation. For flat (closed) background 
geometry for the Universe, the R + R 2 gravity requires a huge initial homogeneous patch 
(huge initial 3-sphere) that contains about 10 11 causally disconnected sub-patches while in 
the R + R 2 supergravity this number is at least 10 6 times smaller. 

We considered topologically trivial FLRW geometries. The homogeneous patch of radius 
Dhomog is enclosed in a smaller volume when k = 1 and in a larger one when k = —1. The 
level of fine tuning is the minimum one when the background spatial geometry has a positive 
curvature. Hence the k = 1 case seems to be favoured unless the Universe is described by a 
non-trivial topology. 

Also, we mention that the study of the pre-inflation supergravitational dynamics revealed 
interesting features such as the initial conditions that give sufficient number of e-foldings, 
that can avoid the eternal process of self-reproduction, and generate a remarkable, however 
ephemeral, anisotropy. 

We underline that throughout this work we focused on pure R + R 2 gravitational and 
supergravitational settings assuming that the dynamics of the higher curvature terms are solely 
responsible for the early Universe inflationary phase. Under these assumptions, our results, 
regarding the initial conditions problem and for a trivial topology, point towards k = 1 and 
theories with extra held dimensions such as the R + R 2 supergravity theory. 
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Figure 10: The figure depicts the number of the causally disconnected regions (# CDR) right after 
the Planck time required to be homogeneous in order for inflation to start. The horizontal axis is the 
time that inflation starts (for the R + R 2 (super)gravity models the Unf is fixed by the CMB with 
tiNF ~ to -1 ~ 10 5 tp\ the plot manifests the initial conditions problem for the low-scale inflationary 
models generally). The solid lines correspond to a flat background spatial geometry, the dotted to 
a positively-curved and the dashed to a negatively-curved one. The red lines correspond to R + R 2 
old-minimal suprgravity, the brown to new-minimal supergravity and the blue to R + R 2 gravity. It 
is shown that the extra supergravity fields decrease the number of CDR. For the k = — 1 case the 
curvature dominates over the fields until the onset of inflation. At H^f = fp the CDR number is 
nonzero; the different initial CDR numbers at Hnf = tp is solely due to the background curvature. 
The CDR number is always larger for the k = — 1 and smaller for the k = 1 background geometry 
regardless the time inflation occurs. 
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